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Abstract 

A general and rigorous method to deal with singularities at the origin 
of a polar coordinate system is presented. Its power derives from a clear 
distinction between the radial distance and the radial coordinate variable, 
which makes that all delta-functions and their derivatives are automatically 
generated, and insures that the Gauss theorem is correct for any distribution 
with a finite number of isolated point-like singularities. 

The method is applied to the Coulomb field, and to show the intrinsic 
differences between the dipole and dimonopole fields in classical electrody- 
namics. 

In all cases the method directly leads to the general expressions required 
by the internal consistency of classical electrodynamics. 



1 Introduction 

How to deal in a consistent and systematic manner with singularities arising at 
the origin of a polar coordinates system is a recurring question in the teaching of 
classical electrodynamics (U |2l [3l |4|| as well as in research [IllOlTl. 

The standard solution is pragmatic: it consists of introducing an appropriate 
delta- function whenever the derivative of an expression that is singular at the origin 
is a distribution rather than zero, i.e., when the derivative is zero everywhere except 
at the origin. Typical examples in three dimensional space are the formulas for the 
divergence of the Coulomb field 



V- A = ^5(|f|)=4vr53(f), (1.1) 
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and the Laplacian of the Coulomb potential 

= V- = -47r53(f). (1.2) 

While these formulas are perfectly correct, as can be verified in a number of ways, 
it would be desirable to have a method such that the "delta-functions" appear at 
the right place without having to remember when this or that expression yields a 
distribution instead of zero. Moreover, one would like to have a technique that is 
general and mathematically rigorous, and whose application is straightforward. 

Indeed, the only thing distribution theory says about an expression that is zero 
everywhere except at a single point derives from the following theorem: A distribu- 
tion which has its support only in one point, say the origin, is a linear combination 
of the delta-function and its derivatives up to a certain order% Therefore, while 
this theorem gives a rigorous justification to the standard method consisting of in- 
troducing the delta-functions "by hand," and of indirectly finding their coefficients 
by a way or another, it does not provide a rule specifying how to formulate an 
initial expression as a distribution so that there is no ambiguity in differentiating 
it in order to get the correct result in a straightforward manner. 

In this paper we present such a method. It is directly applicable to the elemen- 
tary case of an isolated point-like singularity of a scalar or vector field over R^, 
which therefore may be considered as being located at the origin of a spherical 
coordinate system. The method can immediately be extended to n-dimensional 
spaces and fields containing a finite number of isolated point singularities. But it 
cannot easily be extend to topologically more complicated singularities such as, 
for example, the Dirac or Schwinger magnetic monopoles which have a potential 
with a line-like singularity. For similar reasons the distributions considered in this 
paper will always be continuous in the angular variables. 

The essence of the method consists of writing the field under consideration 
in such a way that there is no ambiguity with regards to how to calculate its 
derivatives at every point, including the origin where it may be singular. This 
is made possible by a clear distinction between the radial distance and the radial 
coordinate variable, a distinction that was first made by Tangherlini in the context 
of special and general relativity lUl p. 511-513], which led him to use the sign- 
function sgn(x) to specify how to differentiate at the origin. This idea was 
rediscovered by the author, and others who used it mainly as an ansatz 
whereas it leads in fact to a rigorous method. It is explained and justified in Sec. [21 
where an effort is made to be as correct as possible, but without pretending to a 
rigor such that professional mathematicians would be fully satisfied. In particular, 

'For a proof of that theorem, see |8 , p. 784] or ^ p. 443]. 
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as little mathematical background as possible is used, and the interested reader 
is referred to the numerous publications in which the theory of distributions is 
presented H 

In Sec.[3]it will be shown that Gauss's theorem formulated in polar coordinates 
is correct even if there is a singularity at the origin, provided the method developed 
in this paper is used to calculate the derivatives. This is an important consistency 
check since Gauss's theorem enables to move between the differential and the 
integral formulations of electrodynamics. 

In Sec. |4] the method is applied to the most simple non-trivial case: the 
monopole singularity of a Coulomb or Newton point charge. It will be found that 
an additional singular field appears along side the regular field when it is derived 
from the potential. This additional contribution is essential to insure the local 
conservation of the electromagnetic charge-current density. 

In Secs.|5]and[6]the potentials, fields, and sources of the dipole and dimonopole 
singularities are studied. It will be shown that despite that the dimonopole is 
a scalar-potential-singularity, and the dipole a vector-potential- singularity, their 
fields are identical, except at the origin where the two singularities have different 
delta-function-like fields. It will also be seen that while the derivation of this 
important result is generally somewhat indirect in standard textbooks, it derives 
rigorously from a straightforward application of the present method. 

2 Distributions in spherical coordinates 

In electrodynamics and other areas of physics one is often led to calculating 
integrals of the form 



where T(x) is a well behaved function (i.e., indefinitely differentiable and vanish- 
ing outside a bounded region) and D(x) a function which may be singular, that 
is infinite, at some point. For example T(x) could be the velocity of a stream of 
charges, D{x) an electromagnetic field, and the integral the electromagnetic force 
on that stream. In the language of distribution theory, T is called a test function, 
D a distribution, and the inner product {{D\T)) the value of D on T. 

^For a comprehensive introduction to distribution theory see, e.g., [S] p. 766-798] or ID 
p. 423-541], and for a concise modern presentation [lOJ . 




(2.1) 
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If the distribution is singular at a point xq it is natural to take this point as the 
origin of a polar coordinate system x — xq = r{r, 9, cp). However, unless one is 
very careful, the use of such a coordinate system can lead to ambiguities. For 
instance, the standard practice is to rewrite Eq. (12.11) as follows 

{{D\T))= d(f) define I drr^D{r,e,(f))T{r,e,(l)), (2.2) 
Jo Jo Jo 

where the radial coordinate 

r = \ f\ := Vf-f, (2.3) 

is a positive number. This leads to two difficulties. First, as r G and the 
singularity is at r = 0, that is at the origin of a half-space, a number of basic results 
of distribution theory require some adaptation since they are generally formulated 
for distributions defined over some variable x G M and for singularities at x = 0. 
Second, again as r ^ 0, one has to carefully distinguish between the coordinate 
"r," whose differential dr is part of the integration element, and its magnitude 
"|r|," which may appear in the definition of some particular distribution — which 
for clarity should therefore be written D(r, 9,(1), \r\). 

To better see why this distinction is necessary, let us rewrite Eq. (|2.2I) using 
another spherical coordinate system, parametrized in such a way that the radial 
coordinate is now a signed real number x, i.e., 

{{D\T))= d(f) dOsine dxx'^D{x,e,(l),\x\)T{x,e,(j)). (2.4) 

Jo Jo J ~oo 

In this case the whole set of standard distribution theory theorems can directly be 
used, and there is less temptation to confuse the coordinate x with its absolute 
value \x\. Indeed, while |x| is continuous at 0, its derivative is discontinuous at 
this point since it goes from +1 to —1, or vice versa, when x goes through 0. 
For this reason the second derivative of |s| gives rise to a (5-function at the origin. 
The absolute value |x| is therefore not differentiable at the origin, which is why it 
must be interpreted as a distribution, whereas the coordinate a; is a regular variable 
for which there is no problem at the origin. In fact, the same is true for |r| and 
r in the standard parametrization (12.21) . except that in this case these differences 
are somewhat hidden since the integration over the negative part of the x axis is 
replaced by an extension of the integration range from [0, tt] to [0, 27r]. 

Consequently, to work with distributions in polar coordinates it is necessary to 
insure that all references to a distance or magnitude |x| are made explicit. With the 
spherical parametrization (12.41) . this is best done by introducing the sign function 
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sgn(x) whose definition and relation to Dirac's 5-function are 



sgn(x) = {O X = 0, and sgn(x) = 26(x). (2.5) 

ax 



For any occurrence of |x| we can then write 

d , , , , d\ , 

\x\ = xsgmx), -r\^\ = sgn(x), ^— = 2o(a;), etc., (2.6) 
ax dx^ 

where the rule x5{x) = was used to go from |a;| to d\x\/dx. 

However, the polar parametrization (12.21) is far more frequently used than the 
non-standard parametrization (|2.4I) . In this case we may still use the sign function 
sgn(x) and write |r| = sgn(r) even though r never takes a negative value. But 
this leads to practical difficulties since the radial integration will no more be from 
— oo to +00, but from to +oo, so that, for example, the normalization of Dirac's 
5-function would have to be 1/2 instead of 1 to compensate for the factor 2 in 
Eq. (12.51) . Since this may lead to confusion, and basically consists of formally 
"splitting" the 5-function at r = 0, one may just as well "split" the sign function 
and define a generalized function T(r) such that 



undefined r < 0, 

T(r) := {0 r = 0, and ^T(r) = 6(r). (2.7) 

dr 

r > 0, 



This function should not be confused with Heaviside's step function H(r), which 
is for r < and undefined at r = 0, whereas T(r) is undefined for r < 0, and 
equal to at r = to be equivalent to sgn(x) at x = 0. 



Therefore, for any occurrence of |r| in the standard parametrization (|2.2I) we 
shall write 

d d^ 
\r\ = rT(r), -r\^\ = '^(^)5 = (2-8) 

dr dr^ 

In particular, for any occurrence of the radius vector r we shall write 

f=rT{r)u{9,(p), (2.9) 
where u is the unit vector in the direction of r . 
Finally, we shall make use of the property 

POO 

drT{r)T{r)= / drT{r), (2.10) 
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whose consistency with Eq. (12.71) can be verified by integrating by part the left-hand 
side, i.e., 

X(r)T(-^)(r) -/ dr 5{r)T^-^\r) = drT{r), (2.11) 
Jo Jo 

where we made use of the specification T(0) = 0, as well as of the equation 

dr6{r)T{r) = T{0). (2.12) 



With these rules everything related to the use of a polar coordinate system is 
rigorously taken care of. What remains to make sure is that the generalized func- 
tions D{r, 9, (p, \r\) are really distributions, i.e., that their inner products {{D\T)) 
are converging for any test function T(r) G C°°. This is the case of the distribu- 
tions considered in this paper, which have singularities of type r~" with n > 1 at 
r = 0. Indeed, we must carefully distinguish between the classical functions r^" 
and their derivatives — nr^"^\ which are defined only for r > 0, and the corre- 
sponding distributions which are defined for all r > 0. This is done by defining 
these functions as limits of sequencies of distributions ^Qi p. 51], i.e., 

— := lim— H(r - e). (2.13) 

e— >0 r*^ 



3 Gauss's theorem in polar coordinates 

In order to illustrate the power and generality of the method just explained, we give 
in this section an elementary proof that Gauss's theorem is true in the distributional 
sense for any weakly convergent scalar- or vector-valued distribution F{r) which 
may have a singularity at r = 0. That is, for any finite simply-connected 3-volume 
n bounded by the 2-surface S = dQ, we prove that 



d^j:F{r)= d^nVF{r). (3.1) 




on n 

Referring to Eq. (|2.1I) we have here T = 1, and D = VF{r) is weakly convergent 
because it is a combinations of partial derivatives of the distribution F(r). 

Since F(r) is supposed to be finite and differentiable everywhere except at 
r = 0, the theorem is true by any standard proof for any volume which does 
not contain the origin. This enables to proceed as in any standard analysis of 
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point-like singularities, namely to surround the origin by a small sphere and prove 
the theorem for the leading singularity. However, for simplicity, we postpone to a 
more complete proof the cases where the origin is at the surface of the volume f2, 
and the one where f2 is multiply connected. 

From now one the volume is therefore an ordinary 3-ball of radius R, Oil 
is its surface, and the function F{r) is the leading singularity. If F(r) is simply a 
scalar function /(r), which does not depend on the angular variables, the theorem 
reduces after angular integration to the identity = 0. This remains so after 
multiplying /(r) by a constant vector, so that to deal with non trivial singularities 
we must consider expressions containing the vector r. The most simple such a 
singularity, expressed according to Eq. (12.91) . is 



m = ^ = ^Jir), (3.2) 



which is unbounded at r = for n > 1. We can now operate with V on T(r), and 
make use of Eq. (|2.10l) to remove T(r) where it has become un-necessary, so that 
Eq. (13.11) becomes 



For the 3-ball, we have 

d^i: = d(t) de sin 9 r^u and d^Q = d(j) d9 sin 9 r^dr. (3.4) 
Thus, after an elementary calculation, the three terms appearing in (|3.3I) become 

d^E F{r) = ^irr^u^ 



^ = AnR^-'', (3.5) 




an 

R 



d^n V^—^ =An dr{3- ny^^-""^ = Anr^-" 
Jo 



(3.6) 
(3.7) 

where the lower limits of the volume integrals have not been taken in order to show 




a. 



d^n -^S(r) =471 [ dr r'^-^5{r) = inr^-"" 
Jo ' 



the impact of having used Eq. (12.91) . Indeed, by comparing (|3.6I) and (13.71) . we 
see that the second term on the right hand side of (13.31) has the remarkable effect 
of removing the divergent term at r — > 0, so that the result is equal to the surface 
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integral (13.51) . and, moreover, in the case n = 3, of correcting equation (13.61) so that 
instead of zero it gives the same value as (13.51) . In other words, that second term 
"repairs" the usual formulation of Gauss's theorem in such a way that it becomes 
true for all singularities of the type (13.21) . From there on it is easy to generalize the 
proof to more complicated functions and formulations, so that Gauss's theorem in 
its various forms is true in polar coordinates, even if there are singularities at the 
origin, provided the rules given in Eqs. (|2.7|) to (I2.12|) are followed. 



4 Monopole singularity 

When considering the 1 /r potential of a Newton or Coulomb field, the straight- 
forward application of Eq. (|2.8I) is ambiguous because of the algebraic identity 



(4.1) 



However, in order to obtain the correct form of the field, for which the application 
of Eq. (12.91) is non- ambiguous, it is evident that the potential of a point-charge is 

ipm{^ :=eiT(r), (4.2) 
r 

so that the Coulomb field is 

^™(r) = -V^m = e^T(r) - e^5(r). (4.3) 
The rationalized source-charge distribution is then 

47rp^(f) =V-Em = e^5(r), (4.4) 
which upon integration yields the charge of the source 

q = ipm\l)) = Jjjd'n Pm{^ = e- (4.5) 

There is however an important difference between Eq. (|4.3I) and the usual 
textbook expressions for the Coulomb field: the additional (5-like term on the right 
of (|4.3I) . This requires a careful analysis because the physical interpretation of this 
term is different depending on the physical significance given to the evaluation of 
Em on a test function, i.e., to the inner product {{Em\T)), and to the expression 
Emir) itself. For instance: 
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1. When calculating {{Ern\T)) it is immediately seen that the 5-like term does 
not contribute to the radial integral because the factor in d^il leads to 
a rS{r) product which is zero. Therefore, from the perspective of distri- 
bution theory, the usual Coulomb field Ec = ef/r^ and the field Emir) 
are equivalent, that is corresponding to the same distribution. However, the 
inner product {{Em\T)) has no physical meaning because it has no invariant 
significance, contrary to the charge, Eq. (|4.4I) . The difference between Ec 
and Em has therefore to be evaluated from a physical point of view. 

2. The 5-like term in Eq. (14.31) is a solution of the homogeneous equation 
V ■ Em{r) = 0, which means that it can always be added to a solution of the 
inhomogeneous equation (14.41) . But this does not imply that this singular 
term can always be discarded. For example, in the relativistic case where the 
potential and fields of an arbitrarily moving charge are considered, i.e., in 
the Lienard-Wiechert case, the additional 5-like field is absolutely necessary 
to insure local charge conservation [|71. 

In summary, the derivation of the Coulomb charge distributions (|4.4I) from the 
potential (14.21) has lead to the usual charge distribution, i.e., Eqs. (11.11) or (11.21) . 
whereas the Coulomb field distribution turned out to have the more general form 
(|4.3I) required by the internal consistency of classical electrodynamics. 



5 Dipole singularity 

The dipole singularity, which through extensive experimental verification is found 
to very precisely characterize the intrinsic magnetic dipole moment of elementary 
particles such as the electron, is given by the vector potential 

■= ^'^W' (5.1) 

where \jl \ has the dimension of a charge times a length. The calculation of the 
magnetic field strength gives 

H,{f) = V X I, = (S^l/z-r-) - ^)T(r) + ^^^i^^5(r). (5.2) 

The first term in this expression is well-known, but the one with a (5-function 
is rarely mentioned in textbooks. However, when integrated over 3-space, this 
second term gives the contribution [| TT] p. 184] 

,.s,H.)^^^ = ^/:. (5.3) 
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which is essential in calculating the hyperfine splitting of atomic states fT2l. Thus, 
as in Eq. (|4.3I) for the monopole singularity, we have in Eq. (15.21) a physically 
essential 5-like contribution, which was obtained by a straightforward application 
of the method explained in Sec. |2l 

We can now calculate the sources. As expected, the magnetic charge density 
is zero 

47rp,(f) = V • i?,(f) = 0, (5.4) 
while the rationalized current density is 

47r/,(f) = V X H,{r) = 3^6{r). (5.5) 

Using this current density we can now calculate the magnetic moment by means 
of the standard expression [1 1 , p. 181] to get 

m = ^ /// ^ ^ ^'^^^ ^ 

Therefore, although there are actually no "circulating currents" in the point-like 
distribution (15.51) . the magnetic moment calculated with the formula derived for a 
localized current distribution gives the correct answer. 



6 Dimonopole singularity 

The dimonopole singularity corresponds to the field produced by two magnetic (or 
electric) monopoles of opposite charge ±q separated by an infinitesimal distance 
|A|. The potential for such a field is therefore the scalar expression 

Vdmir) ■= - - — (6.1) 



r 



|f + A| 



At large distance, or at vanishingly small separation A, we can take for this potential 
the first term of the Taylor development, i.e., 

fdm{r) ^ g^(A ■ f)T(r). (6.2) 

The field strength is then 

HU^ = -V^,rn = (34(/^ ■ ^1 - 4)^W - ^^^^W, (6.3) 
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where we have defined 

fI=qX. (6.4) 

Expression (16.31) is remarkably similar to the corresponding expression (15.21) for 
an intrinsic dipole, and it can be seen that the difference between a dipole and a 
dimonopole field is entirely contained in the point-like singularity at the origin, 
i.e., 

H,{rl - HUrl = ^-^^^^5{r) + ^-^6{r) = ^5(r). (6.5) 

As a result, when integrated over 3-space, the dimonopolar 5-singular term in 
Eq. (16.31) gives the contribution |[TT1 p. 141] 

-///d'n3^i(r) = -f/:, (6.6) 

which differs in sign and in magnitude from the corresponding expression (15.31) 
for an intrinsic dipole. It is this difference which enables to conclude that the 
dipolar fields from distant stars are produced by magnetic dipoles, rather than by 
magnetic dimonopoles [fT2ll . 

We can now calculate the sources. As expected, the current density is zero 

^7ijdm{r) = V X Hdm{r) = 0, (6.7) 
while the rationalized charge density is 

47rp,^(rO = V • Hdmi^ = S^S{r), (6.8) 

i.e., a distribution that is odd in r so that the total charge is zero, as it should be for 
a dimonopole. We can finally calculate the first moment of this charge density by 
means of the standard expression for a charge distribution [fTTl p. 137]. This gives 




(6.9) 



a result formally similar to Eq. (15.61) which was obtained by using the formula for 
the moment of a current distribution. This illustrates again that despite the great 
similarity of their fields at a distance from the origin, the dipole and dimonopole 
singularities are in fact very different. 
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